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The choice for parametric techniques in the dis- 
cussion article is motivated by the claim that for 
multivariate extreme- value distributions, "owing to 
the curse of dimensionality, nonparametric estima- 
tion has essentially been confined to the bivariate 
case" (Section 2.3). Thanks to recent developments, 
this is no longer true if data take the form of multi- 
variate maxima, as is the case in the article. A wide 
range of nonparametric, rank-based estimators and 
tests are nowadays available for extreme-value cop- 
ulas. Since max-stable processes have extreme- value 
copulas, these methods are applicable for inference 
on max-stable processes too. The aim of this note 
is to make the link between extreme-value copulas 
and max-stable processes explicit and to review the 
existing nonparametric inference methods. 

1. EXTREME-VALUE COPULAS 

Let the random variables Y\ , . . . , Yd represent the 
maxima in a given year of a spatial process (e.g., 
rainfall) that is observed at a finite number of sites, 
x±,...,xd, m a region X in space W (typically, 
p = 2). Let Fi,...,Fd be the marginal cumulative 
distribution functions, assumed to be continuous. In 
the article, these are assumed to be univariate gen- 
eralized extreme-value distributions, an assumption 
that will not be needed here. 

The random variables Ud = FdiYd) are uniformly 
distributed on the interval (0, 1) and the joint cumu- 
lative distribution function C of the vector U\,. . . , 
Ud is the copula of the random vector Y\ , . . . , Yd : 

(1) C(«i, . ..,u D ) = Pr(C/i <u 1 ,...,U D < u D ) 
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for < Ud < 1. The requirement that the random 
vector Y\, . . . , Yd is max-stable entails 



(2) C m (u 1 1 /m ,...,u D i m ) = C(u 1 ,...,UD) 

for all m > 0. In [18], it was shown that (2) holds if, 
and only if, 

(3) C(ui,...,u D ) = exp{-rA(vt, . . . ,v D )}, 



l/m> 



_1 log Ud- The 



,UD) 

where r = — Yld=i ^°& u d and Vd - 
domain of the Pickands dependence function A is 
the unit simplex, Sd = {v G [0, l] D : Yld v d = !}■ 
A necessary and sufficient condition for a function A 
on (Srj to be a Pickands dependence function is that 

A{vx, ...,Vd) 

(4) 

= / max(uisi,...,UDSD)dM(si,...,SD) 
Js D 

for a Borel measure M on Sd verifying the con- 
straints f s s<jdM(si, . . . , sd) = 1 for all d € {1, ... , 
D}. In particular, A is convex and max(ui, . . . , vd) < 
A(vi, . . . , vd) < v\ + • • • + vd- In dimension D = 2, 
these two properties completely characterize Pickands 
dependence functions (but not if D > 3). 

2. MAX-STABLE MODELS 

The representation in (3)-(4) is valid for general 
max-stable copulas and therefore also holds for the 
finite-dimensional distributions of the max-stable 
processes considered in Section 6 in the article. The 
purpose of this section is to make this relation ex- 
plicit. 

Consider the simple max-stable process 
(5) Z(x) = max[Sj max{0, Wj(x)}], x G W, 

where {<S 3 }|2 =1 are the points of a Poisson process 
on with rate s~ 2 ds and where W\,W2, ■ ■ ■ are 
iid replicates of a stationary stochastic process W 
on MP, independent of the previous Poisson pro- 
cess, and such that E[W + (x)] = 1, where we write 
W + (x) = max{0, W(a;)}. Particular cases of this 
model include the so-called Smith model [24], the 
Schlather model [22] and the Brown-Resnick mod- 
el [12]. 
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The stationary, marginal distribution of Z(x) in (5) 
is unit-Frechet and the joint distribution function of 
the vector Z(xi), . . . , Z(xd) is given by 

Pt[Z(x{) <zi,..., Z(x D ) < z D ] 



exp 



-/'■ 



({(s,w) : max(w + (x d )/z d ) > 1/sj 



for z d > 0, where [i is the intensity measure of the 
Poisson point process {(Sj,Wj)}'jL 1 . A simple cal- 
culation shows that 

Pt[Z(xi) <zi,..., Z{x D ) < z D ] 

= exp( — E max{W + (x d )/z d } 
\ L d 

As a consequence, the copula of Z{x\), . . . , Z(xd) 
is given by the extreme- value copula with Pickands 
dependence function 



(6) A(v u 



,VD) 



E 



max{v d W + (x d )} 



for (vi, . . . ,vd) G Sd- As illustrated by the compu- 
tations in [5], the integral arising in (6) can rarely 
be calculated analytically, unless D = 2. 

To recover the spectral measure M in (4) from 
the distribution of the stochastic process W, let R = 
^2 d W + (x d ). On the event R > 0, consider the ran- 
dom vector {W + {xi),...,W + {x D ))/R. Then 

dM{ Sl ,...,s D ) 

= Pr(R > 0)E[R | Vd: W + {x d ) /R = s d ;R> 0] 

• Pr [Vd : W + (x d ) /R G ds d \ R > 0] . 

3. ESTIMATION 

Nonparametric estimators of the Pickands depen- 
dence function are surprisingly easy to construct and 
calculate. The starting point is the simple fact that 
for u € [0, 1] and for (v±, . . . , vd) € So, the extreme- 
value copula C with Pickands dependence function A 
satisfies 

(7) C(u Vl ,...,u VD )=u A ( vl '-< VD \ 

as can be verified from (3). Using (7), the function A 
can be recovered from the copula C in various ways, 
for instance, through integrals of the form 



(8) 



f(C(u v \...,u VD ))g(u)du 



f(u a )g(u)du, 



for well-chosen functions / and g. Plugging estima- 
tors for C and solving for a then yields estimators 
for A. 

A natural estimator for C is the empirical copula. 
Let (Yn, . . . , Yin), for i € {1, ... , n}, be an indepen- 
dent random sample from a distribution with con- 
tinuous margins and copula C. The empirical copula 
is defined as 



C n (ui, . . . ,u D ) 



(9) 



1 n 

- V/{F n i(y<i) < «i, . . . ,F nD (Y w ) < u D }, 
n ' 



i=l 



a = A(v 1 ,...,v D ) 



where F nd is the (marginal) empirical distribution 
function of Yi d , . . . , Y nd . Being based on multivariate 
ranks, the empirical copula is invariant under mono- 
tone transformations of the data. The empirical cop- 
ula goes back to the seminal paper by Riischendorf 
[21] and has been studied and applied intensively, 
such as recently in [23, 25, 26]. 

Inserting the empirical copula into (8) and solv- 
ing for A(vi, . . . ,Vd) produces simple and (almost) 
explicit estimators. Particular instances are the 
Pickands estimator [6, 8, 11, 19] and the Caperaa- 
Fougeres-Genest estimator [4, 8, 11]. The bivariate 
versions of these estimators are special cases of the 
weighted estimator in [17]. Minimum-distance esti- 
mators are another instance of this technique [3]. 
Standard errors can be obtained via resampling [2, 
20] or via empirical likelihood [17]. 

A drawback of the nonparametric estimators of A 
is that they typically do not produce valid Pickands 
dependence functions — remember the representation 
in (4) that such functions must satisfy. A way to 
overcome this issue is by projecting a possibly in- 
valid pilot estimator A n onto the family of Pickands 
dependence functions, yielding 

(10) Ar j = argmin/ (A n - A) 2 , 

AeA D Js D 

where Ad denotes the family of all Pickands de- 
pendence functions in dimension D and where the 
integral is with respect to some measure on So- In 
general, the infinite-dimensional least-squares prob- 
lem in (10) does not admit an explicit solution. Ap- 
proximate solutions can be obtained by perform- 
ing the minimization over the (finite-dimensional) 
class of Pickands dependence functions with discrete 
spectral measures supported on a given, finite grid 
[7, 11]. 
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Finally, note that a nonparametric estimator A n 
can be transformed into a parametric one via mini- 
mum-distance or projection techniques: in (10), re- 
place Ad by the parametric model of interest. If the 
model happens to be specified via the point process 
representation (5), then this technique requires the 
calculation of the Pickands dependence function A 
via (6). 

4. TESTING 

Nonparametric methods are particularly suitable 
for hypothesis testing. Of special interest are the hy- 
pothesis of max-stability in general and the good- 
ness of fit of a parametric model in particular. In 
most cases, critical values are computed via resam- 
pling methods. 

Even if the data at hand are vectors of component- 
wise maxima, it is a good idea to test whether it 
is safe to assume that the underlying distribution 
is max-stable, in particular, when the end-goal is to 
perform prediction and/or extrapolation. For bivari- 
ate extreme- value copulas, the first two moments 
of the random variable W = C(U U U 2 ) = F(Y U Y 2 ) 
happen to satisfy a particular linear relation. The 
sample moments of the random variables W n \ , ■ • • , 
W nn defined by 

n 

W ni = - V I(Y tl < Y a ,Y t2 < Y i2 ) 
n 

t=i 

can therefore be converted to a test statistic for the 
null hypothesis of max-stability [1, 10]. 

Another approach for testing max-stability is by 
comparing the empirical copula C n in (9) with the 
extreme-value copula that has a given estimator A n 
as its Pickands dependence function. For the bi- 
variate case, Cramer-von Mises tests based on the 
Pickands and Caper aa-Fougeres-Genest estimators 
are described in [3, 14]. 

Finally, the adequacy of the hypothesis of max- 
stability can be tested by directly exploiting the cop- 
ula max-stability relation (2) through a comparison 

of C n (u\, . . . ,ud) with C™(u^ m , . . . , u 1 ^ 171 ) for var- 
ious values of m > 0. Cramer-von Mises type test 
statistics turn out to be particularly effective [16]. 

The goodness of fit of a parametric model can 
be tested by comparing the fitted parametric es- 
timator for A with a nonparametric one [9]. For 
max-stable models arising through the point pro- 
cess representation in (5), the function A has to 
be computed through the relation (6). Shape con- 
straints such as exchangeability can be tested simi- 
larly [15]. 



5. CONCLUSION 

Nonparametric methods yield an attractive alter- 
native inference method for max-stable dependence. 
Estimators and test statistics of the Pickands de- 
pendence function are (almost) explicit, even in the 
general, multivariate case. Moreover, as the proce- 
dures are based upon the ranks of the data only, the 
step of modeling the margins can be skipped (which 
is not to be confused with the false statement that 
the uncertainty on the margins has been eliminated 
altogether). Many of the methods described in this 
contribution are implemented in the R package cop- 
ula [13]. 
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